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1.  INTRODUCTION 


One  of  the  important  parameters  in  shell  design  is  the  total  aerodynamic  drag.  The  total  drag  consists 
of  three  components:  the  pressure  drag  or  the  wave  drag  (excluding  die  base),  die  viscous  drag,  and  the 
base  drag.  The  base  drag  component  is  a  large  part  of  die  total  drag  and  can  be  as  high  as  50%  or  more 
of  die  total  drag.  Of  these  three  components  of  drag,  die  most  difficult  one  to  predict  is  die  base  drag 
because  it  depends  on  the  pressure  acting  on  die  base.  Therefore,  it  is  necessary  to  predict  die  base 
pressure  as  accurately  as  possible. 

The  ability  to  compute  the  base  region  flow  field  for  projectile  configurations  using  Navier-Stokes 
computational  techniques  has  been  developed  over  the  past  few  years  (Sahu,  Nietubicz,  and  Steger  1985; 
Sahu  1986,  1987).  Recently,  improved  numerical  predictions  (Sahu  and  Steger  1988;  Sahu  1990;  Sahu 
and  Nietubicz  1990)  have  been  obtained  using  the  Cray-2  supercomputer  and  a  more  advanced  zonal 
upwind  flux-split  algorithm.  This  zonal  scheme  preserves  die  base  comer  and  allows  better  modeling  of 
die  base  region  flow.  These  studies  have  included  base  flows  for  different  base  geometries.  This 
capability  is  very  important  for  determining  aerodynamic  coefficient  data,  including  the  total  aerodynamic 
drag.  As  indicated  earlier,  a  number  of  base  flow  calculations  have  been  made,  and  base  drag  and  total 
drag  have  been  predicted  with  reasonable  accuracy.  However,  because  available  data  are  lacking,  die 
predictive  apabilities  have  not  been  assessed  with  detailed  base  pressure  distributions,  mean  flow  velocity 
components,  and  turbulence  quantities.  This  is  especially  true  of  base  flow  for  axisymmetrical  bodies  at 
transonic  and  supersonic  speeds.  Recently,  experimental  measurements  (Herrin  and  Dutton  1991)  have 
been  made  in  die  base  region  for  supersonic  flow  over  a  cylindrical  afterbody.  The  data  include  base 
pressure  distribution  (along  the  base),  mean  flow,  and  turbulence  quantities. 

Figure  1  is  a  schematic  diagram  showing  the  important  features  of  supersonic  base  flow.  The 
approaching  supersonic  turbulent  boundary  layer  separates  at  die  base  comer,  and  the  flee  shear  layer 
region  is  formed  in  die  wake.  The  flow  expands  at  the  base  comer  and  is  followed  by  die  rccompression 
shock  downstream  from  the  base  that  realigns  the  flow.  The  flow  then  redevelops  in  the  trailing  wake. 
A  low  pressure  region  is  formed  immediately  downstream  from  the  base,  which  is  characterized  by  a  low 
speed  recirculating  flow  region.  Interaction  between  this  recirculating  region  and  the  inviscid  external 
flow  occurs  through  the  flee  shear  mixing  region.  This  is  die  region  where  turbulence  plays  an  important 
role. 
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Expansion  Waves 


Figure  1.  Schematic  Diagram  of  Supersonic  Base  Flow. 

The  basic  configuration  used  in  this  study  is  a  cylindrical  afteibody.  As  mentioned  earlier,  a  simple 
composite  grid  scheme  has  been  used  for  accurate  modeling  of  the  base  comer.  Numerical  flow  field 
computations  have  been  performed  at  M_=  2.46  and  at  (f  angle  of  attack.  Three  turbulence  models  (two 
algebraic  models  and  a  two-equation  model)  are  used  in  the  base  flow  region.  All  the  amputations  have 
been  performed  on  the  Cray-XMP  supercomputer.  Details  of  the  flow  field  such  as  Mach  number 
contours  and  base  pressure  distributions  are  presented.  Computed  base  pressure  distributions  are  compared 
with  available  experimental  data  for  the  same  conditions  and  the  same  configuration.  The  algebraic 
turbulence  models  predict  a  large  change  in  the  base  pressure  distribution  over  the  base.  The  two-equation 
k-e  model  predicts  a  rather  small  change  in  base  pressure  along  the  base  and  compares  very  well  with  the 
experimentally  measured  base  pressure  distribution. 

2.  GOVERNING  EQUATIONS  AND  SOLUTION  TECHNIQUE 

The  complete  set  of  time -dependent,  Reynolds-averaged,  thin  layer  Navier-Stokes  equations  is  solved 
numerically  to  obtain  a  solution  to  this  problem.  The  numerical  technique  used  is  an  implicit,  finite 
difference  scheme.  Although  time-dependent  calculations  are  made,  the  transient  flow  is  not  of  primary 
interest  at  the  present  time.  The  steady  flow,  which  is  the  desired  result,  is  obtained  in  a  time  asymptotic 
fashion. 


2.1  Governing  Equations.  The  complete  set  of  three-dimensional  (3-D),  tune-dependent,  generalized 
geometry.  Reynolds-avenged,  thin  layer,  Navier-Stokes  equations  for  general  spatial  coordinates  r\,  and 
C  can  be  written  as  (Pulliam  and  Steger  1982) 

♦  dn<j  *  d;/?  -  Re‘'di$  (1) 

in  which 

4  *  4C*.  y.  *.  /)  -  longitudinal  coordinate;  q  =  tj(jc,  y,  z.  t)  -  circumferential  coordinate 
(  =  C(x,  y,  z,  t )  -  neatly  normal  coordinate;  t  =  t-  time 


pU 

PuU  +  %xp 
pvU  + 
pwU  *  \,p 
{e+p)U  -  l,p 


pV  pW 

puV  +  x\j>  puW  *  t>jp 

C  -  i  pvV  +  Ti^  #«i  pvW+C/r 

J  J 

pwV  ♦  tl,p  pwW  +  t>j> 

(e+p)V-i\tp  (e*p)W-^,p 
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and  in  which 


i-i 

/ 


*  £  *  C,2>»{  *  *  t,»{  * 

MC2  *  ^  ♦  £>*{  *  *<(,«<  *  c,v;  .  C,wc)C, 

♦  ?5  ♦  £>n  *  *  Vc  *  C“t>t, 


<?i  *  £  ♦  c,2) 


JL(«2  +  v2  +  w2)^ 


Kfl; 


/»Ky-UJ 


+  ~(Cx«  ♦  CyV  +  C,w)(^«c  +  +  c,wc) 


(3) 


In  Equation  1,  the  thin  layer  approximation  is  used  and  the  viscous  terms  involving  velocity  gradients 
in  both  the  longitudinal  and  circumferential  directions  are  neglected.  The  viscous  terms  are  retained  in 
the  normal  direction,  £,  and  are  collected  into  die  vector  §.  These  viscous  terms  are  used  everywhere. 
However,  in  the  wake  or  the  base  region,  similar  viscous  terms  are  also  added  in  die  stream-wise  direction. 
For  this  computation,  the  diffusion  coefficients  p  and  K  contain  molecular  and  turbulent  parts.  The 
turbulent  contributions  are  supplied  through  either  algebraic  or  a  two-equation  k-e  turbulence  model. 

The  velocities  in  the  Ij,  t],  and  £  coordinate  directions  can  be  written 

t/  =  ^+  u^+v^  +  w^ 

Vsilj  +  iiTk+vtly  +  wn, 

W^St  +  u^  +  v^  +  w^ 

which  represent  the  contravariant  velocity  components. 
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The  Cartesian  velocity  component!  (a,  v,  w)  are  retained  as  the  dependent  variables  and  are 
nonuimensionalized  with  respect  to  a^  (die  free  stream  speed  of  sound).  The  local  pressure  is  determined 
using  the  relation 

P  ■  (T  -  1  > [ «  -  0.5p(u2  v2  w2)]  (4) 

in  which  y  is  the  ratio  of  specific  heats.  Density,  p,  is  referenced  to  pM  and  the  total  energy,  e,  to  pM  a^2. 
The  transport  coefficients  are  also  nondimensionalized  with  respect  to  the  corresponding  free  stream 
variables.  Thus,  die  Prandtl  number  that  appears  in  Sis  defined  as  Pr  =  cfjtjKm.  In  differencing  these 
equations,  it  is  often  advantageous  to  difference  about  a  known  base  solution  denoted  by  subscript  0  as 

-  4>)  +  V*  '  4>)  +  8n(^  -  60)  *  8;(rf  -  fi0)  -  Re~lS^S  -  S0) 

*  -3t(5o  “  ~  ^5^0  +  R*  (5) 


in  which  6  indicates  a  general  difference  operator,  and  d  is  the  differential  operator.  If  die  base  state  can 
be  properly  chosen,  the  differenced  quantities  can  have  smaller  and  smoother  variation  and  therefore  less 
of  a  differencing  error  (Pulliam  and  Steger  1982). 

2.2  Numerical  Technique.  The  implicit,  approximately  factored  scheme  for  the  thin  layer  Navier- 
Stokes  equations  using  central  differencing  in  the  ti  and  ^  directions  and  upwinding  in  £  is  written  in  die 
following  form 


[/  +  hfyA'f  +  h%C*  -  hRe'1Z(iJ-1tiHJ  -  £>;|;] 

x  [/  ♦  h^(A~)*  *  hSn6*  -  D i  |nJA0" 

-  -Ar(8f [(/*♦)"  -O  +  -C]  +  -O 

♦  -  O  -  Re  -1  8C($ "  -  Sm))  -De(6H-6m)  (6) 
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in  which  h  -  Ar  or  (Ary 2  and  the  free  stream  base  solution  is  used.  Here,  6  is  typically  a  three  point 
second  order  accurate  central  difference  operator,  $  is  a  midpoint  operator  used  with  the  viscous  terms, 
and  the  operators  5^b  and  are  backward  and  forward  three-point  difference  operators.  The  flux  F  has 

AAA  A 

been  eigensplit  and  the  matrices  A,  B,  C,  and  M  result  from  local  linearization  of  the  fluxes  about  the 
previous  time  level.  Here  /  denotes  the  Jacobian  of  the  coordinate  transformation.  Dissipation  operators 
De  and  Di  are  used  in  die  central  space  differencing  directions.  The  smoothing  terms  used  in  die  present 
study  are  of  the  form: 

D,\n  -  ♦  e46^i^.S3J|ny 

-  (Ar)/“1[e25p(5)p8  +  2.5e45p(il)5]|t,y 


in  which 

ft-  Ifr»l 

Id  I 

and  p(B)  is  die  true  spectral  radius  of  B.  The  idea  here  is  that  the  fourth  difference  will  be  tuned  near 
shocks  (e.g.,  as  P  gets  large,  the  weight  on  the  fourth  difference  drops  down  while  the  second  difference 
tunes  up). 

For  simplicity,  most  of  the  boundary  conditions  have  been  unposed  explicidy  (Sahu  1987).  An 
adiabatic  wall  boundary  condition  is  used  on  the  body  surface  and  the  no-slip  boundary  condition  is  used 
at  the  wall.  The  pressure  at  the  wall  is  calculated  by  solving  a  combined  momentum  equation.  Free 
stream  boundary  conditions  are  used  at  the  in-flow  boundary  as  well  as  at  the  outer  boundary.  A 
symmetry  boundary  condition  is  imposed  at  the  circumferential  edges  of  die  grid  while  a  simple 
extrapolation  is  used  at  die  downstream  boundary.  A  combination  of  symmetry  and  extrapolation 
boundary  condition  is  used  at  the  center  line  (axis).  Since  the  free  stream  flow  is  supersonic,  a 
nonreflection  boundary  condition  is  used  at  the  outer  boundary.  The  flow  field  is  initially  set  to  free 
stream  conditions  everywhere  and  then  advance  in  time  until  a  steady  state  solution  is  obtained. 

2.3  Composite  Grid  Scheme.  In  the  present  work,  a  simple  composite  grid  scheme  (Sahu  1990)  has 
been  used  in  which  a  large  single  grid  is  split  into  a  number  of  smaller  grids  so  that  computations  can  be 
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performed  on  each  grid  separately.  These  grids  use  the  available  cm  memory  one  grid  at  a  time.  The 
remaining  grids  are  stored  on  an  external  disk  storage  device  such  as  the  solid  state  disk  device  (SSD)  of 
the  Cray  X-MP/48  computer.  The  Cray-2  has  a  large  in-core  memory  to  fit  the  large  single  grid. 
However,  for  accurate  geometric  modeling  of  complex  projectile  configurations,  which  include  Mum  noses, 
sharp  comets,  and  base  cavities,  it  is  also  desirable  to  split  the  large  data  base  into  a  few  smaller  zones 
on  the  Cray-2  as  well.  The  use  of  a  composite  grid  scheme  requires  special  care  in  storing  and  fetching 
die  interface  boundary  data  (Le.,  the  communication  among  die  various  zones).  In  die  present  scheme, 
there  is  a  one  to  one  mapping  of  die  grid  points  at  die  interface  boundaries.  Thus,  no  interpolations  are 
required.  Details  of  die  data  storage,  data  transfer,  and  other  pertinent  information  such  as  metric  mi 
differencing  accuracy  at  the  interfeces  are  given  in  the  work  of  Sahu  and  Steger  (1987)  and  Sahu  (1988). 

2.4  Turbulence  Modeling.  For  the  base  flow  calculations,  three  turbulence  models  have  been  used. 
Two  of  these  are  algebraic  eddy  viscosity  models  (Baldwin-Lomax  model  and  Chow  model).  The  third 
one  is  a  two-equation  k-e  turbulence  model  which  is  also  an  eddy  viscosity  model 

2.4.1  Baldwin-Lomax  Model  This  model  is  die  one  developed  by  Baldwin  and  Lomax  (1978).  It 
is  a  two-layer  model  in  which  an  eddy  viscosity  is  calculated  for  an  inner  and  an  outer  region.  The  inner 
region  follows  the  Prandtl-Van  Driest  formulation,  hi  both  die  tamer  and  outer  formulations,  the 
distribution  of  vorticity  is  used  to  determine  the  length  scales,  thereby  avoiding  the  necessity  of  finding 
the  outer  edge  of  die  boundary  layer.  For  die  tamer  region, 

<P<)«*r-P*2l®l  (7) 

in  which 

I  -  iey[l-exp(-y  VO] 
y*m(Pw« ■  /(VPw> 


and  I  to  I  is  the  absolute  magnitude  of  vorticity.  The  eddy  viscosity  for  die  outer  region  is  given  by 


O*|)on*r  "  ^^epP^wokt^Utb^y^ 


(8) 
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in  which  FlljJU  *  ymax  Fnn,  or  ymax  u^/Fmrr.  the  smaller  of  the  two  values.  The  quantities 
and  Fmtrr  are  detennined  from  the  function  F(y)  *  ylcol  [1  -  exp  (-y*  /  >4*)],  in  which  is  the 
maximum  value  of  F(y)  and  is  die  value  of  y  at  which  it  occurs.  The  function  F,J»  is  the 
Klebanoff  inteimittency  factor.  The  quantity  is  the  difference  between  die  maximum  and  minimum 
total  velocity  in  die  profile  and,  for  boundary  layers,  the  minimum  is  zero. 

The  outer  formulation  can  be  used  in  wakes  as  well  as  in  attached  and  separated  boundary  layers.  Bor 
free  shear  layer  flow  regions  or  wakes,  die  Van  Driest  damping  term  [exp{-y*  /  V)]  is  neglected.  Also, 
for  the  base  or  wake  region,  the  distance  y  is  measured  from  the  center  line  of  symmetry.  It  is  necessary 
to  specify  the  following  constants:  A*  =  26,  Ccp  =  1.6,  =  0.3,  =  0.25,  K  =  0.4,  and  K  -  0.0168. 

This  type  of  simple  model  is  generally  inadequate  for  complex  flows  containing  flow  separation  regions 
such  as  base  flow. 

2.4.2  Chow  Model.  Another  algebraic  model  that  has  been  used  in  some  of  our  base  flow 
computations  is  that  of  Chow  (1985).  This  model  is  intended  to  be  used  in  the  base  or  wake  region  only. 
It  is  based  on  the  simple  exchange-coefficient  concept  The  turbulent  eddy  viscosity  coefficient  is  usually 
given  by 

*  "  7^2  x“*  (9) 

in  which  x  is  die  distance  measured  from  the  origin  of  the  mixing  region  (i-e.,  the  base).  ue  is  the  velocity 
at  the  edge  of  the  mixing  region,  and  o  is  the  spread  rate  parameter.  It  is  known  that  a  assumes  a  value 
of  12  for  incompressible  flow  and  it  increases  slightly  with  Mach  number. 

O  =  12  +  2.76  Me 

in  which  Me  is  given  by 


8 


be  as  ceitained  from 


As  a  first  approximation,  the  average  value  of  p,  is  assumed  to  be  same  at  all  points  for  a  constant  x 
location.  After  leattachment,  turbulence  should  decay.  Since  the  interest  in  the  base  flow  calculations 
is  to  obtain  the  correct  base  pressure,  it  is  assumed  that  the  eddy  viscosity  level  at  the  reattachment  nays 
the  same  at  other  locations  downstream.  For  base  flow  with  jet,  similar  algebraic  relation  %  used 
for  the  jet  shear  layer.  This  model  as  well  as  Baldwin-Lomax  (1978)  model  are  algebraic  models  and 
depend  only  on  local  information.  The  two-equation  model  contains  less  empiricism  and  allows  die  flow 
history  to  be  taken  into  account. 

2.4.3  Two-Eauation  k-e  Model  The  two-equation  turbulence  model  used  here  is  Chien’s  (1982)  k-e 
model  which  is  similar  to  that  of  Jones  and  l  aunder  (1972).  In  tins  model,  two  transport  equations  are 
solved  for  the  two  variables,  k  (turbulent  kinetic  energy)  and  e  (turbulent  dissipation  rate). 


Dk  afm  ]a*  [  a**,-  d«,- 

pof  *  ^  k  *  * 


-  Pe  "  2p~ 

y* 
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Here,  yn  is  the  distance  normal  to  the  surface.  The  coefficients  in  the  k  and  e  equations  are  given  by 

Cj  -  1.44  c3  -  10.  ok  •  1.0,  oe  -  1.3 

c2  •  1.921 1  -0.3exp(-J?/)]  c f  -  0.09[1  -exp( -0.01y*)] 

in  which  Rt  *  i^/vc. 


The  1-c  model  employs  the  eddy  viscosity  concept  and  relates  die  turbulent  eddy  viscosity  to  k  and  e  by, 

P,  -  cpp(*2/e)  .  (12) 

Following  the  same  procedure  used  for  the  mean  flow  equations,  die  turbulence  field  equations  can  be 
written  in  conservation  form  and  then  transformed  into  generalized  coordinates  (Sahu  and  Danberg  1986). 
The  resulting  axisymmctric  set  of  transformed  turbulence  equations  cm  be  written  as 


3d.  di.  3(5.  i 
3t  3£  3d  Re 


3/?,  dft  . 


(13) 


in  which 
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and  P  is  the  production  ten  given  m 

r  -  Mi  *i*  £)(«;  +  vc  ♦  w?>  *  p»(C«“c  +  Cyvc +  Z*w02 

*  Pt(£  *  i  *  £)(“*  +  ^  +  w£)  +  M%x^  ♦  ♦  ^H^)2 

This  constitutes  a  low  Reynolds  number  formulation  of  the  k-e  model  Calculations  are  extended  to  the 
wall  itself,  and  exact  values  of  the  dependent  variables  at  the  wall  are  used  as  boundary  conditions. 
Oden's  model  is  better  mathematically  behaved  near  the  wall  and  is,  dais,  used  in  this  study. 


The  turbulence  field  equations  are  solved  using  the  Beam  and  Warming  (1978),  implicit, 
approximately  factored,  finite  difference  scheme.  A  convenient  solution  algorithm  for  these  equations  has 
die  following  sequence: 

[(f  -  AID  ")  +  Al(8^B  "  -  6^C  ")]  A$,"  -  RHS  ( 13)  (14a) 

[/  ♦  Ai(5^A ■  -  %Nm)]  Aq*  -  Af,"  (14b) 

«,"*1  -  *,*  .  A (Me) 

in  which  RHS  is  die  right-hand  ride  of  Equation  13.  A  and  B  are  Jacobian  matrices  resulting  from  the 
local  linearization  of  die  flux  terms  Et  and  Gv  The  source  terms  are  treated  implicitly.  This  results  in 
the  Jacobian  matrices  C,  N,  and  D,  which  are  included  in  the  £  and  £  operators  as  shown  in  Equation  14. 
The  Jacobian  matrices  are  given  as 
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in  which 


(c*.  *  ti  * 


The  (pernors  and  85  are  central  difference  operators.  The  numerical  smoothing  is  based  on  an  up¬ 
wind  scheme,  and  die  details  are  given  in  Sahu  (1984). 

3.  MODEL  GEOMETRY  AND  EXPERIMENT 

The  computational  accuracy  of  a  numerical  scheme  can  be  established  through  comparisons  with 
available  experimental  data.  The  model  used  in  the  experiment  and  in  the  computational  study  is  shown 
in  Figure  2.  It  is  an  axisymmetrical  cylindrical  afterbody,  which  has  a  diameter  of  63.5  mm.  This  figure 
also  shows  the  stations  where  mean  and  fluctuating  velocity  components  were  measured  with  a  Laser 
Doppler  Velodmcter  (LOV)  system.  The  same  configuration  is  used  in  die  numerical  simulations  for  a 
direct  comparison. 


Experimental  measurements  (Herrin  and  Dutton  1991)  for  this  model  have  been  made  at  the  University 
of  Illinois  supersonic  wind  tunnel.  The  model  was  tested  at  0°  angle  of  attack,  Mach  number  of  2.46,  and 
Reynolds  number  of  521  x  107  per  meter.  In  addition  to  measuring  the  velocity  components  at  a  few 
selected  longitudinal  positions  in  the  wake  or  base  region,  the  base  pressure  was  measured  at  19  positions 
along  the  base.  Such  detailed  base  pressure  measurements  have  not  been  made  in  the  past  and  are  very 
helpful  in  the  code  validation  process.  The  velocity  profile  is  also  measured  at  a  station  upstream  from 
the  base,  which  provides  the  upstream  boundary  condition  for  base  region  flow  field  calculations. 

4.  RESULTS 

Numerical  computations  have  been  made  for  d *  cylindrical  afterbody  at  a  Mach  number  2.46  and  at 
0°  angle  of  attack.  The  three-plane  version  of  the  3D  code  was  run  for  die  0°  angle  of  attack  case.  Two 
end  planes  woe  used  to  specify  symmetrical  boundary  conditions  in  the  circumferential  plane. 

The  solution  technique  requires  the  discretization  of  die  entire  flow  region  of  interest  into  a  suitable 
computational  grid.  The  grid  outer  boundary  has  been  placed  1  diameter  away  from  die  surface  of  the 
afterbody.  The  downstream  boundary  was  placed  10  diameters  away  from  the  base.  Since  the  calculations 
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Figure  2.  Afterbody  Measurement  Locations. 

ire  in  the  supersonic  regime,  the  computationil  outer  boundary  was  placed  close  to  the  body  and  a  no 
reflection  boundary  condition  used  at  that  boundary.  Figure  3  shows  an  expanded  view  of  the  grid  in  die 
bare  region.  The  surface  points  on  die  afterbody  and  the  base  were  obtained  first  These  were  then  used 
as  inputs  for  obtaining  die  full  grid  using  an  algebraic  grid  generation  program.  The  full  grid  is  split  into 
two  zones,  one  upstream  from  the  base,  and  the  other  one  in  the  base  region  or  the  wake.  These  grids 
consist  of 22  x  60  and  95  x  119  grid  points,  respectively.  Figure  3  shows  the  longitudinal  grid  clustering 
near  the  base  coiner.  Grid  points  are  also  clustered  near  the  afletbody  surface  to  capture  the  viscous 
effects  in  the  turbulent  boundary  layer.  These  clustered  grid  points  are  spread  out  downstream  of  the  base 
in  the  wake  to  capture  the  free  shear  layer  region.  For  the  0°  angle  of  attack  case  considered,  the  grid  was 
rotated  circumferentially  5°  on  either  side  of  the  midplane.  This  provided  the  three  planes  needed  in  the 
code  to  use  central  finite  differences  in  the  circumferential  direction.  In  each  case,  the  solution  was 
marched  from  free  stream  conditions  everywhere  until  the  final  converged  solution  was  obtained.  The 
results  are  now  presetted  for  both  mean  and  turbulence  quantities.  Comparison  of  the  computed  results 
is  made  with  the  available  experimental  data  (Herrin  and  Dutton  1991). 
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A  few  qualitative  results  are  presented  next  Figure  4  shows  the  pressure  contc«r  plot  for  the  base 
region.  The  features  to  observe  are  the  flow  expansion  at  the  base  comer  followed  by  the  recompression 
shock  downstream  from  the  base  (coalescence  of  contour  lines).  Figures  Sa  and  5b  show  die  comparison 
of  the  computed  Mach  number  contours  with  experimentally  obtained  Schliren  photograph  of  the  base 
region  flow  field.  Both  the  experiment  and  the  computed  results  show  the  flow  expansion  at  the  base  and 
the  recompression  shock  downstream  from  the  base.  In  addition.  Figures  5a  and  5b  show  the  free  shear 
layer  in  the  near  wake.  Although  not  indicated  in  Figure  5a,  the  flow  in  the  near  wake  is  primarily 
subsonic.  Figure  6  shows  the  computed  velocity  vectors  in  the  base  region.  The  recirculatory  flow  in  the 
near  wake  is  clearly  evident.  Flow  reattachment  occurs  at  about  three  base  radii  downstream  from  the 
base.  The  magnitude  of  the  velocity  is  shown  to  be  quite  small  in  the  immediate  vicinity  of  die  base.  The 
computed  results  shown  in  Figures  4, 5,  and  6  were  obtained  using  die  two-equation  k-e  model. 

Figures  7  and  8  show  the  velocity  components  in  the  stream-wise  and  normal  directions,  respectively. 
These  velocity  profiles  are  taken  at  four  longitudinal  positions  in  the  wake  or  die  base  region  (X/D  =  1.26, 
1.42, 1.73,  and  1.89).  The  computed  velocity  profiles  obtained  using  two  algebraic  turbulence  models  and 
the  two-equation  k-e  model  are  compared  with  the  experimental  data.  Figure  7  shows  the  comparison  of 
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tiie  u  (stream-wise)  component  of  velocity.  In  general,  the  profiles  obtained  with  the  k-c  model  agree 
much  better  with  the  experiment  in  the  shear  layer  regions  for  X/D  =  1.26  rod  X/D  =  1.42.  The  profiles 
are  rather  poorly  predicted  by  both  the  algebraic  models  at  these  two  stations.  The  reattachment  point 
estimated  from  the  experimental  measurements  is  about  1.4  base  diameters  downstream  from  foe  base. 
The  computed  value  with  the  k-e  model  is  1.5.  This  small  disagreement  is  also  seen  in  the  flow 
redevelopment  region  downstream  from  the  reattachment  (X/D  =  1.73  and  1.89).  The  algebraic  turbulence 
models  predict  the  reattachment  point  better  than  that  predicted  by  the  k-e  model  The  velocity  profiles 
predicted  with  these  models  agree  fairly  well  with  foe  experimentally  obtained  profiles  at  these  two 
stations.  Chow  model  predictions  are  slightly  better  thro  those  predicted  by  the  Baldwin-Lomax  model 
in  this  flow  redevelopment  region.  Figure  8  shows  the  comparison  of  the  w  (vertical)  component  of  the 
velocity.  This  component  of  velocity  is  better  predicted  by  the  k-e  model  than  foe  algebraic  models  both 
in  foe  flow  recirculation  and  redevelopment  regions.  The  profiles  by  foe  algebraic  models  do  not  agree 
well  with  the  experimental  data  for  radial  positions  greater  than  half  of  tire  base  radius. 

Some  of  the  turbulence  quantities  are  presented  next  Figure  9  shows  the  turbulent  kinetic  energy 
profiles  at  the  same  longitudinal  positions  in  the  wake.  The  computed  k  profiles  are  obtained  using  the 
two-equation  k-e  turbulence  model  In  foe  recirculation  region  (X/D  =  1.26)  and  near  the  reattachment 
(X/D  =  1.42),  the  peak  observed  experimentally  (r/R  =  0.5)  is  poorly  predicted  by  the  k-e  model 
(r/R  =  0.4).  The  location  and  foe  magnitude  of  the  peak  agree  somewhat  better  at  X/D  =  1.42  than  at 
X/D  =  1.26.  The  agreement  of  the  computed  profiles  with  the  data  is  good  in  the  flow  redevelopment 
region  (X/D  =  1.73  and  1.89).  Figure  10  shows  the  turbulent  dissipation  rate  (e)  profiles  at  die  same 
positions  in  the  wake.  As  seen  in  this  figure,  e  increases  from  the  center  line  of  symmetry  with  radius 
to  about  0.4  of  the  base  height  where  the  peaks  occur  and  then  drops  quickly  to  very  small  values  at  about 
r/R  =  0.65.  The  magnitude  and  foe  location  of  foe  peaks  decrease  smoothly  with  increasing  axial  distances 
downstream  from  the  base. 

Figure  11  shows  foe  turbulent  shear  stress  profiles  in  the  wake.  The  computed  values  obtained  by 
both  the  algebraic  models  and  the  k-e  model  are  compared  with  the  experimental  data.  In  general,  a  small 
improvement  can  be  observed  in  foe  predicted  values  with  the  k-e  model  over  the  algebraic  models. 
Discrepancy  exists  between  the  experimentally  obtained  turbulent  shear  stress  and  the  predicted  shear 
stresses  with  all  the  turbulence  models.  This  is  true  especially  near  the  peaks  at  X/D  =  126  and  1.42. 
The  magnitude  of  foe  peak  predicted  by  the  k-e  model  is  about  the  same  as  predicted  by  the  Baldwin- 
Lomax  model  at  these  two  positions;  however,  they  both  underpredict  the  experimental  peak.  The  Chow 
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model  underpredicts  the  peak  even  more.  As  lbr  the  location  of  fhe  peak,  the  k-e  model  does  better  fan 
the  algebraic  models.  As  X/D  is  increased  from  1.26  to  1.42,  the  location  of  die  peak  predicted  by  the 
k-e  model  moves  closer  to  die  center  line  similar  to  that  observed  in  the  experiment  This  is  not  seen  in 
the  prediction  by  the  algebraic  models.  The  k-e  model  predictions  agree  better  than  the  predictions  by 
the  algebraic  models  at  X/D  *  1.73  and  1.89. 

Of  particular  interest  is  die  accurate  prediction  or  determination  of  base  pressure  and,  hence,  base  drag. 
Figure  12  shows  die  base  pressure  distribution  (along  the  base).  The  base  pressures  predicted  by  both  the 
algebraic  models  and  the  two-equation  k-e  turbulence  model  are  compared  with  the  experimental  data 
(Herrin  and  Dutton  1991).  The  experimental  data  are  shown  in  dark  circles,  and  the  computed  results  are 
shown  in  lines.  Here,  Z/D  =  0.0  corresponds  to  the  center  line  of  symmetry  and  Z/D  =  05  corresponds 
to  the  base  comer.  The  base  pressures  predicted  by  both  algebraic  turbulence  models  show  a  big  increase 
near  the  center  line  of  symmetry.  The  experimental  data  show  almost  no  change  (only  3%)  in  the  base 
pressure  distribution.  The  base  pressures  are  very  poorly  predicted  by  the  algebraic  models,  not  wily  near 
the  center  line  but  also  near  the  base  craner.  A  much  improved  base  pressure  distribution  is  predicted  by 
the  k-e  model,  and  its  agreement  with  the  measured  base  pressure  is  quite  good.  The  k-e  prediction  shows 
a  small  increase  in  the  base  pressure  near  the  center  line,  which  is  not  observed  in  the  data. 
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Figure  12.  Base  Pressure  Distribution.  M_  «  2.46.  a  =  0. 
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5.  CONCLUDING  REMARKS 


A  zonal,  implicit,  tune-marching  Navier-Stofces  computational  technique  has  been  used  to  compute 
the  turbulent  supersonic  base  flow  over  a  cylindrical  afterbody.  Flow  field  computations  have  been 
performed  at  M,*,  =  2.46  and  at  the  angle  of  attack,  a  =  0.0.  Various  eddy  viscosity  turbulence  models 
(two  algebraic  and  a  two-equation  k-e)  have  been  used  to  provide  the  turbulence  closure.  The  k-e 
equations  were  formulated  in  a  generalized  coordinate  system  and  were  solved  using  an  implicit  algorithm. 

Numerical  results  show  the  details  of  the  flow  field  such  as  Mach  number  contours,  pressure  contours, 
and  velocity  vector  plots.  Comparison  of  both  tire  mean  and  turbulence  quantities  has  been  made  with 
the  available  experimental  data.  Both  algebraic  turbulence  models  predict  the  mean  velocity  components 
poorly  in  the  recirculatory  flow  region  in  the  wake.  In  general,  tire  velocity  components  predicted  by  the 
two-equation  k-e  model  agree  better  with  the  experimental  data  than  tire  algebraic  models  do.  Discrepancy 
exists  between  the  predicted  turbulent  shear  stress  and  the  experiment  for  all  these  turbulence  models. 
A  small  improvement  in  die  predicted  location  and  magnitude  of  tire  peak  in  shear  stress  exists  with  the 
k-e  model.  Computed  base  pressure  distributions  have  been  compared  with  tire  measured  base  pressures. 
The  base  pressures  predicted  by  the  algebraic  models  show  a  much  larger  variation  and  do  not  agree  well 
with  the  data,  compared  to  the  k-e  model.  The  measured  base  pressures  show  a  very  small  change  along 
the  base  and  are  predicted  rather  well  with  the  k-e  turbulence  model. 
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